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1. Introduction 

Our present knowledge of strong interactions indicates that QCD (Quantum Chro- 
modynamics) is the quantum field theory that models the physics of this fundamental 
interaction. However, many important properties of the strongly interacting parti- 
cles, the hadrons, can not be described using only perturbative QCD. This happens 
because the QCD coupling is large at low energies. Thus, in order to calculate static 
properties of hadrons, or to describe their structure, one needs additional tools. 

In the last years, many interesting models to study non perturbative aspects 
of hadronic physics were developed based on the idea of gauge/string duality. The 
main source of inspiration for these kind of models was the discovery of the AdS/CFT 
correspondence [1, 2, 3] that is an exact duality between string theory in certain ten 
dimensional geometries and supersymmetric SU(N) gauge theories with large N on 
the corresponding boundary. In particular, string theory in AdS 5 x S* 5 space is dual 
to a four dimensional gauge theory. 

In the AdS/CFT correspondence, the gauge theory is conformal. The idea of 
breaking this conformal invariance by an Infrared (IR) cut off, represented in the 
dual AdS geometry by a geometrical cut off in the radial coordinate of the space, 
was introduced in [4] as a tool to reproduce the high energy scaling of hadronic 
scattering amplitudes for processes at fixed angles. This scaling had been found in 
QCD a long time before [5, 6] but the corresponding string theory description was 
lacking for a very long time. 

This approach of considering a maximum radial size of AdS space as the dual 
of an infrared cut off in the gauge theory was then used in [7, 8] to calculate the 
mass spectrum of glueballs. In these articles, boundary conditions were imposed on 
fields living in an AdS slice and the corresponding normal modes were associated 
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with hadronic states. This kind of model, later called AdS/QCD hard wall, was then 
applied to other hadrons, as, for example, in [9]. 

The hadronic masses calculated using the hard wall model do not exhibit the 
observed approximate linear relation between mass squared and spin (and radial 
excitation number), corresponding to the so called linear Regge trajectories. Rather, 
the hard wall masses show asymptotically quadratic trajectories. This motivated a 
different AdS/QCD approach consisting of a background involving AdS space plus a 
field that acts effectively as a smooth infrared cut off [10]. This, so called soft wall 
model, leads to linear Regge trajectories for mesons (see also [11] for glueballs). 

For fermions, the soft wall model, as originally formulated [10], does not lead to a 
discrete spectrum because the dilaton introduced in the action is factorized out in the 
equations of motion. This situation can be changed, following alternative approaches 
in order to reproduce Regge behavior for baryons. One example can be found in 
[12, 13], where the authors consider models in asymptotically AdS space including 
a warp factor in the metric. Other possibility, studied in [14, 15, 16], considers a 
z dependent (or dressed) mass for the fermionic modes propagating in AdS space. 
This last approach has been considered in [15] to study nucleon form factors, and 
in [17] to obtain some generalized parton distributions (GPD) for nucleons. Here 
we consider baryons in an AdS/QCD soft wall model with dressed mass in order to 
describe deep inelastic scattering (D1S). 

The study of deep inelastic scattering (DIS) using the AdS/CFT correspondence 
appeared first in [18]. Then other authors considered the description of DIS using 
various AdS/QCD models like, for example, in [19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 
29]. In [19] the soft wall AdS/QCD model was considered for scalar hadrons and a 
hybrid model involving a soft wall cut off for the photons and a hard wall cut off 
for the fermions was also discussed. However, the calculation of baryonic structure 
functions using a soft wall cut off was still missing. We will use the modified soft wall 
model with dressed mass, studied in [15, 16] to calculate the DIS structure functions 
for baryons. Then we compare our F 2 with experimental results for the proton 
in the regime of large values of the Bjorken parameter x where the supergravity 
approximation that we use is more reliable. 

In section 2 we present a briefly review of DIS and hadronic structure functions. 
In section 3 we show our AdS/QCD calculation of the structure functions. Then in 
section 4 we analyze our results, discussing the possible choices of the parameters 
of the model, plotting the structure function F 2 for some kinematical regimes and 
comparing with experimental results. 

2. Deep Inelastic Scattering and Structure Functions 

Deep inelastic scattering (DIS) is a process where a highly energetic lepton, an elec- 
tron in general, interacts with a hadronic target through the exchange of a virtual 



-2- 



p 



X 



Figure 1: Diagram of a deep inelastic scattering. 



photon, as show in the diagram of figure 1. The momenta of the photon and of the ini- 
tial hadron are respectively and P M . After the interaction there is a final hadronic 
state represented by X with momentum P£- The experimental measurement of the 
inclusive cross section of DIS corresponds to detecting the final lepton, thus deter- 
mining the momentum transfer q^, but not the final hadronic state X. That means, 
summing over all possible final hadronic states X. One usually parametrizes DIS 
using as dynamical variables the photon virtuality q 2 and the Bjorken parameter 



Deep inelastic scattering, in the strict sense, corresponds to the limit q 2 — > oo, with 
x fixed. 

The inclusive cross section for DIS can be calculated from the hadronic tensor 
that is defined as 



where J^{y) is the electromagnetic hadronic current and s the spin of the initial 
hadron. For an unpolarized scattering (spin independent) this tensor can be decom- 
posed in terms of the structure functions Fi(x,q 2 ) and F 2 (x,q 2 ) as [30] 



W>" = Fl (x, q 2 ) (rT - \f ) + %F 2 (x, g 2 ) (P" + |) (P* + q -) . (2.3) 



Considering only final states with just one baryon with mass M x , we can intro- 
duce a basis of such one particle states and write the hadronic tensor as 



x = — 



q 2 /2P-q. 



(2.1) 



W " V = / Ae^(P,*| \ny),r(0)] \P,s), 



(2.2) 
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= \ E E 5 { M x + ( P + ?)* ) < P > + 9» ( P + sx|^(0)|P, s) . 

s,sx M x 

(2.4) 

So, in order to calculate the hadronic structure functions one needs to find the 
matrix elements of the current and also the spectrum of hadronic masses M x ■ We 
will present in the next section the calculation of these quantities using the modified 
soft wall model. 



3. DIS in the soft wall model with fermionic dressed bulk 
mass 

In the soft wall model, a four dimensional gauge theory is represented by a gravity 
dual consisting of fields living in anti-de Sitter space AdS$, whose metric can be 
written as 



ds 2 ee g mn dy m dy n = ^(dz 2 + ^dy^dy") = %{dz 2 - dt 2 + {dyf) , (3.1) 

z A z* 

with the presence of an additional background of the form e' 1 ^^ . This factor plays 
the role of an infrared cut off, with the energy parameter k representing the cor- 
responding scale. The general form of an action integral corresponding to a given 
Lagrangian density C is 

I = J dzd 4 y e^ 2 C . (3.2) 

For the case of scalar and vector fields, this kind of action integral leads to 
dual four dimensional theories with infrared cut off and discrete mass spectra. For 
fermions one must do some modification in order to find a dual theory with an 
infrared cut off. This happens because using a standard fermionic Lagrangian in the 

2 2 

action above the background e~ K z factors out of the equations of motion. So, the 
fermions are not affected by the energy scale k and get no discrete mass spectrum. 
This problem can be solved introducing a z dependent mass, as studied in [15, 16]. 
Following this approach, the appropriate action that describes the dynamics of the 
fermionic and gauge fields and their interaction is 



/ 



dzd 4 y \f—g e K 



' TP TP r - 



if 

2 



(3.3) 
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where the field strength is F mn = d m A n — d n A m and e™ = <5™-| and the covariant 
derivative is: 



with: = ~\{^ h m — 5 b z 5m) ■ The coupling g v will be associated with the electric 
charge e and T a = (7^, — 27 s ) are the Dirac matrices. 

In order to produce a discrete spectrum with the Regge behavior for the fermions 
we introduced the effective fermionic potential Vp{z) = k 2 z 2 /R depending on the z 
coordinate (associated with the boundary energy scale), to dress the fermionic bulk 
mass. 

The physical motivation for this procedure is the following. On one hand, accord- 
ing to the AdS/CFT dictionary, the mass of a supergravity bulk field is related to the 
dimension of the corresponding boundary operator. On the other hand, in general 
the dimensions of quantum operators receive anomalous contributions, depending on 
the energy scale. Since the energy scale of the boundary theory is holographycally 
related to the localization in the z coordinate, the possibility of anomalous contribu- 
tions to the dimension of the operators can be translated into z dependent masses for 
the dual bulk modes [14, 31]. This kind of procedure makes it possible to introduce 
an important ingredient of QCD that is not considered in most of the AdS / QCD 
models. Other related works that consider masses varying in the bulk can be found 
in, for example, [32, 33, 34, 35]. 

For the gauge fields it is convenient [18] to impose the gauge condition 



Note that from now on we use the notation of raising and lowering the four 
dimensional indices with the Minkowski metric: A^ = rf v ' A v and □ = ri^d^dy . 
With the gauge choice (3.5) the equations of motion that emerge from the action 
(3.3) are 



We impose the condition that the boundary value of the gauge field represents a 
virtual photon with polarization rf and space-like momentum 



V m = d m --u^[T a ,T b }- igv A 



(3.4) 




(3.5) 




(3.6) 



.m-v 



(3.7) 



The corresponding solutions are 



A,(z, y) = rip e*™ k 2 r(l + f^) z 2 U{\ + 
A z (z,y)= l -r ] -qe^ T(l + ^)zU(l + 




(3.8) 
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where U(a; b; w) are the confluent hypergeometric functions of the second kind. 

Now, regarding the fermionic fields, it is convenient to do the following field 
redefinition 



y(y,z) = e + * 2z2 ^(y,z). 
in such a way that the equations of motion take the form 



+ 7 5 <9 2 



1 
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i 2 2 

m + k z 



(3.9) 



(3.10) 



where jd 



dfj_. The quantity /i, the bulk fermion mass was replaced by the 



dimensionless parameter m = fiR. In order to solve the equation of motion, we 
decompose the fermionic field into chiral components 



il>(y,z) = ip L (y,z)+ip R (y,z) , ij} L/R 

with ( 7 5 ) 2 = 1. 

We will consider solutions of the form 



1 TT 



(3.11) 



1>l/r(v,z) = e ip y l -(l Tl b )u s {P) ^h/ R {z) , 



(3.12) 



that contain the plane wave factor and the four component spinor u s {P) correspond- 
ing to a four dimensional free fermion with momentum P M and spin s. Then, the z 
dependent parts of chiral components of the fields: /l/r(^) must satisfy 



-dl + «V + 2k" (rn T 1 -) + m( ^ 2 ±1} ] f L/R (z) = -P 2 f L/R (z) . (3.13) 

Following the usual prescription of gauge/gravity dualities, normalizable solu- 
tions for fields in the gravity side are dual to states in the boundary four dimensional 
theory. Equations (3.13) have normalizable solutions only when — P 2 , the four di- 
mensional mass, has the discrete values 

- Pi = M 2 n = Ak 2 (u + m + ^) , (3.14) 

with n — 0, 1, 2, The corresponding discrete set of normalizable solutions f2/ R ( z ) 
are 



f2( Z ) = / 2r (^+l) m + 3/2 m+1 £ -^/2 L m+l/2 (/c 2 2) (3^5) 

w T(n + m + 3/2) v ; ' K J 



f»( Z ) = / 2r (^+l) m+1/2 m -*V/2 L m-l/2 ( 2 2) (g^g) 

,/fiV ; y T(n + m + 1/2) « v ; v ; 
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with normalization condition f dz f™, R (z)f2/ R (z) = S mn . 

o 

We want to describe processes where the initial state is a proton that absorbs a 
virtual photon transforming into a final state corresponding to an excited hadronic 
state of spin 1/2. So we consider the interaction action 



S mt [h X}=g v J dzd^y^e'^ 2 | A m * x 7 m % . (3.17) 

where represents the initial proton, that we take as the state with lowest mass 
level, corresponding to n — 0. So the initial momentum Pi = p satisfies — p 2 = Mq = 
4k 2 (m + 1/2). The fermionic field represents a final state with (higher) mass 
Mx and momentum Px = p + q satisfying —P x = M\ = 4k, 2 (nx + to + 1/2), where 
we are representing as nx the integer associated with the excitation level of the final 
state. The corresponding solutions have the form 



(3.18) 

where Sj and sx are the spins of the initial and final fermionic states. 

In order to find out the structure functions for the hadrons we have to calculate 
the interaction action (3.17) with the solutions for the fields. We can simplify the 
calculations by considering, as it was done in [18, 19], that we are probing the hadron 
with a particular photon with polarization n^ satisfying n ■ q = 0. For this situation 
the z component of the gauge field does not contribute and, substituting the solutions 
(3.6) and (3.18) in (3.17) the interaction action, for this initial and final states (i,X) 
takes the on shell form 

S int [i,X] = ^ (2nY5\Px-p-q)v\n Sx (P x )^{-^)u s Xp)Tdn x ) 

+ u Sx (P x )r(^Y^)u Si (p)^RM] , (3.19) 



where T L (n x ) and T R (n x ) are integrals involving the two fermionic solutions with 
different chiralities. They have a similar structure and can be written, in terms of 
the variable w = k 2 z 2 , in the general form: 



2 poo 2 

l(m,n x ) = C(m,n x )T(l + f- 2 ) J dww^ 1 e- w U(l + 2; w) , (3.20) 
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where 

y 1 (m — 1)1 (nx + m — 1) 

The integrals X^(n x ) and Tn(n x ) correspond to Z(m, nx) with m = m + 5/2 and 
rh = m + 3/2 respectively. Performing the integral (3.20) we get 

^rM /<~r" r(jl + , ix ) 

Z(m,n x ) = ^- v (3.22) 

2/t 2 r(m-l)r^ + nx + mJ 

Now that we calculated the action in the soft wall background for photons inter- 
acting with fermions with dressed mass, we connect this result with the four dimen- 
sional boundary theory using a proposal similar to the one used in refs. [18, 19]. In 
these references the matrix element of the fermionic electromagnetic current on the 
boundary four dimensional theory was considered to be equal to the corresponding 
bulk interaction action. Here we will take a different point of view and consider that 
bulk/boundary duality implies that these quantities are proportional, rather than 
necessarily equal. So we assume the relations: 



V ,{P x \J^q)\P) ={2Tif5\P x -P-q)^{P + q\.P{Q)\P) = JC eff S int [h 

^(P\J»(q)\P x ) = (2tt) 4 5\P x -P-q)r ] ,(P\r(0)\P + q) = K.„ S int [X, i] , 

(3.23) 

where )C eff plays the role of a bulk/boundary effective factor that phenomeno logically 
adjust the bulk supergravity quantities to the boundary observed ones. Following this 
prescription, the hadronic tensor in eq. (2.4), contracted with the photon polarization 
i], can then be written in terms of our interaction action of eq. (3.19) as 



v,v,w^ = \ ml + (p + g) 2 ) % x 



4 ^ v * v , 4 

M x 



Si s x 

+l 2 L u Sx rr ( - ) u s M Si rr { - ) u Sx + X 2 R u Sxl ^ + \ Si u Si YT^u Sx } (3.24) 

where we defined g e ff = fC ef} g v and = and omitted the dependence of 
the spinors on the momenta. Note that we are summing over the final spins and 
averaging over the initial ones. Using the property 

^2(u a ) a (p)(u 8 ) p (p) = (7% + m) aP , 
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satisfied by both the initial and final state (on shell) spinors, with the corresponding 
masses and momenta, we find 



V.VuW^ = \ J2 S ( M x + (P + <l) 2 ) 9e ff { -X 2 L (n x )X 2 R (n x )M x M oV - V 

M x 

+ (x 2 L (n x ) + X 2 R {n x )) ((p • V? ~ \(V 2 + P ■ q) V ' v) } (3-25) 

The sum over the final states X can be approximated by an integral over a 
continuum of states as it was done in in refs. [18, 19]. In the present case we this 
corresponds to replacing the sum over X of the delta functions by the factor: 
That means: 



VtfuW'™ » ^{-X 2 L (n x )X 2 R (n x )M x M rj ■ rj 

l 2 (n x ) + X 2 R (n x )) ((p • v) 2 ~ \(P 2 + P ■ q) V • ^ } (3.26) 



For the particular photon that we are considering, with 77 • q = we get from eq. 



(2.3) 



2r 1 2 4(VP) 2 F 2 



(3.27) 



Comparing this with our expression for the hadronic tensor eq. (3.26) we find 
our results for the fermionic structure functions 



Fi 



9eff 

16k 2 



X L (n x ) +X R (n 



2 4x 



-X L {n x )X R {n x )M \ M 2 + 



q 2 (l -x) 



x 



9eff 

32k 2 



{x 2 L (n x )+X 2 R (n x )) 



Si 

x 



(3.28) 



The excitation level n x of the final state is not an independent variable. It can 
be expressed in terms the DIS variables q 2 and x. Using eqs. (2.1) and (3.14) one 
finds 



-iP + q) 2 = M 2 + q 2 



1 



- 1 



x 



so 



n x = 



An 2 



M 2 X =AK 2 (n x + m + ^j, (3.29) 

(3.30) 
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q 2 [GeV 2 ] 

Figure 2: F 2 as a function of q 2 , where g 2 is in units of (GeV) 2 . The plot consider 
x = 0.85. The higher curve in for m = 0.6, the middle one for m = 0.7 and the lower is 
plotted with m = 0.8. Dots correspond to experimental data [36, 37]. 

4. Analysis of the results 

In order to calculate the structure functions and analyze their dependence on x 
and q 2 , first we have to fix the parameters of the model. The effective coupling 
g e fj = K eff gv contains the electric charge g v , that satisfies g 2 = 1/137 but mul- 
tiplied by the yet undetermined parameter K,,. inserted in the model in order to 
phenomenologically adjust the relation between bulk and boundary quantities and 
essentially fix the size the structure functions. We found that the choice g 2 ^ = 0.66 
leads to structure functions F 2 with values compatible with experimental data. So 
we used this value in the calculations. 

The infrared energy scale k is related with the Regge slope. In fact, in the present 
holographic model the Regge slope is 4k 2 . We will choose 4k 2 = 0.9 GeV 2 that was 
found in ref. [12, 14] to give a good adjustment for the trajectory of nucleons. So 
we use k = 0A74GeV. The other parameter of the model: m, is considered as a free 
parameter, associated, as we discussed before, with the anomalous dimension of the 
operator that creates the baryonic states. We used this free parameter to look for 
the best fit of the shape of the structure function F 2 to the experimental data. It is 
important to remark that the supergravity approximation used in the model is, in 
principle, valid only for large values of x. This happens because, as discussed in [18], 
for low values of x string theory corrections would become relevant. So we analyzed 
the region 0.8 < x < 1.0. We show in Fig. 2 the structure functions found with the 
holographic model for x = 0.85 compared with the corresponding experimental data 
that appear in PDG [36, 37] considering the values of m that led to the best fits, 
that means, m close to 0.7. 

In order to check if the dependence of F 2 on q 2 has a similar form for other 
values of x in the range considered, we show in Fig. 3 this structure function for 
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Figure 3: F 2 as a function of q 2 . The plot on the left consider x = 0.8 and on the right 
x = 0.9. In both plots the higher curve in for m = 0.6, the middle one for m = 0.7 and the 
lower is plotted with m = 0.8. 

x = 0.8 and x = 0.9. In both cases we included plots for three different choices of 
the parameter m. We see that F 2 decreases with q 2 in a way that is compatible with 
the experimental results shown in figure (16.7) of [36] or in [37] for x = 0.85. 

Now, let's consider the dependence of F 2 on the Bjorken parameter x. There are 
not many experimental results for this structure function for x close to one. But one 
finds in ref. [38] an interesting investigation of F 2 at large x, for low values of q 2 , 
using some parametrizations obtained from experimental data. In particular, they 
show results for the F 2 as a function of x for the cases of q 2 = 4 GeV 2 and 9 GeV 2 . 
Their results show a decrease in this structure function when we increase x in the 
interval 0.9 < x < 1. 

We show in Fig. 4 the structure function F 2 obtained using our holographic 
model with z dependent mass in the range 0.8 < x < 1.0, for these two values of q 2 
analyzed in [38], using values for m around 0.7. The order of magnitude of our results 
was adjusted by the choice of the effective coupling g e ff to be consistent with the 
experimental values, so they are of the same order of those found in [38], however, 
in contrast to the results of [38], we found an increase in the structure functions 
when x — > 1. So, our model for DIS of baryons does not give a good description 
for the dependence of F 2 on x at low q 2 . This may be a consequence of the fact 
that in this simple model the final hadronic states have just one baryon with spin 
1/2. In a general non elastic process the final state may include more than one 
hadron and also hadrons with higher spins. A description of such kind of process 
with multiple hadrons and also higher spins is out of the scope of the present kind 
of model. Nevertheless it is interesting to see that the model is able to reproduce, 
by adjusting the free parameters, the experimental dependence of F 2 on q 2 . 

Finally, it is interesting to consider the high q limit of the structure functions 
eq.(3.28). Using the property 
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Figure 4: F 2 v/s x. The plot at the left consider q 2 = 4GeV 2 and the other is for 
q 2 = 8GeV 2 . In both plots the higher curve in for m = 0.6, the middle one for m = 0.7 
and the lower is plotted with m = 0.8. 



T(a + y) / r l^ 6 ~ B ^ l /1( -j 



we can expand the integrals I^m denned by eq. (3.22) and that appear in eq.(3.28). 
This way we find properties shared with other AdS/QCD models at dominant order 
in q, like 

F 2 = 2F X (4.2) 
that implies that for x — > 1 we approach the Callan-Gross relation. We also find 



V4k 2 / 



x m+ i(l-x) m+ K (4.3) 



g z = 30 GeV ' 
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Figure 5: F 2 v/s x. The plot at the left consider q 2 = 30 GeV 2 and the other is for 
q 2 = 1000 GeV 2 . In both plots the higher curve in for m = 0.6, the middle one for m = 0.7 
and the lower is plotted with m = 0.8. It can be seen that at high q the structure functions 
decrease for x — > 1. 
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Notice that when m = r — 3/2 our model reproduces results found in previous 
works that consider DIS in holographic models [18, 19]. Additionally, at high values 
for q 2 , we find a decreasing behavior for the structure function when x is close to 
one, as can be seen in Fig. 5. 
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